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Abstract: Results from the multiplicity distributions of inclusive photons and charged 
particles, scaling of particle multiplicities, event-by-event multiplicity fluctuations, and charged- 
neutral fluctuations in 158-A GeV Pb-|-Pb collisions are presented and discussed. A scaling 
of charged particle multiplicity as Np^J^^^'^^ and photons as Np^^^^'^^ have been observed, 
indicating violation of naive wounded nucleon model. The analysis of localized charged-neutral 
fluctuation indicates a model-independent demonstration of non-statistical fluctuations in both 
charged particles and photons in limited azimuthal regions. However, no correlated charged- 
neutral fluctuations are observed. 

1. Introduction 

The year 2000 has been quite interesting and remarkable for heavy ion physicists search- 
ing for a new form of matter, the quark-gluon plasma (QGP). In the beginning of the 
year, CERN has announced that there is definite evidence of possible QGP signals from 
the experiments performed at the SPS [0]. This declaration was made after scrutinizing 
results from all the dedicated experiments which have been taking data since 1994 using 
the Pb beams at 158 GeV/nucleon (a total energy of 33 TeV). In the middle of the year, 
this field entered a new era with the commissioning of the Relativistic Heavy Ion Collider 
(RHIC) at BNL which is a captive machine for QGP search. Soon after, new results from 
the collisions at RHIC have started to appear ||^. At this time, it is most appropriate 
to review one of the basic observables of the nuclear collisions, that is the particle mul- 
tiplicity distributions. One of the distinct advantages of the heavy ion collisions at such 
high energies is the production of large number of particles in every event. This allows 
for a detailed study of event-by-event fluctuations in particle multiplicities and ratios of 
identified particles. This is important for understanding the evolution of nuclear system 
at high energy collisions. Formation of disoriented chiral condesates (DCC), which is a 
direct consequence of chiral symmetry restoration, would also give rise to event-by-event 
correlated fluctuations in charged particles to neutrals. In this manuscript we discuss 
these topics in view of the data from the WA98 experiment at the CERN-SPS. 

For a thermalized system undergoing a phase transition, the variation of the tempera- 
ture with entropy density is interesting as the temperature is expected to increase below 
the transition, remain constant during the transition, and then increase again P,^. This 
can be studied by measuring the mean transverse momentum, (pr), and particle multiplic- 
ities for varying centrality, for a number of colliding systems at different energies. Since 
multiplicity of produced particles is an important quantity to characterize the evolving 
system, its fluctuation from event to event may provide a distinct signal of the phase 
transition from hadronic gas to QGP phase. Recently, several new methods using event- 
by-event fluctuations of (px) and particle multiplicities have been proposed to probe into 
the mechanism of phase transition |^|-^. Another interesting phenomenon is the forma- 
tion of domains of DCC p|-p!2| which gives rise to isospin fluctuation, in which the neutral 
pion fraction can deviate significantly from 1/3, the value for uncorrelated emission of 
pions. Several methods have been proposed to search for signals of DCC [|T^-[T^. The 
most direct signal comes from the event-by-event fluctuations in the number of charged 
particles to photons in localized (rj-cp) phase space. 

The observed fluctuations will have contributions from statistical fluctuations and those 
which have dynamical origin. The contribution from dynamical origin comprise of (a) 
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Figure 1. Minimum bias multiplicity distributions of charged particle and photon-like 
clusters for Pb+Pb reactions at 158- A GeV. The other three curves show multiplicity 
distributions for three different centralities, 0-1%, 0-2% and 0-5%. 

fluctuations which do not change event-by-event, e.g., those from Bose-Einstein (BE) 
correlations, resonance decays and (b) the fluctuations which have new physics origin 
and may vary from event-to-event. The major interest for us would be to probe into the 
fluctuations which have new physics origin, such as those arising near the tricritical point 
of QCD phase diagram and the formation of DCC To extract the dynamical part from the 
observed fluctuations, one has to understand the contributions from statistical and other 
known sources. It is possible to probe at the non-statistical fluctuations from experimental 
data in a model independent way by comparing these with mixed events generated from 
the data. Once properly understood, mixed events provide the best means to infer about 
the presence of non-statistical fluctuations. 

This manuscript is organized as follows. In the next section, we discuss multiplicities of 
charged particles and photons and their pseudorapidity distributions. In section 3 we dis- 
cuss the scaling of particle multiplicities with centrality of the reaction expressed in terms 
of the number of participants. In section 4 we present multiplicity fluctuations and discuss 
the importance of making proper centrality selection. In section 5 we discuss charged- 
neutral fluctuations in the full phase space of the multiplicity detectors. In section 6 we 
present the results of an analysis to search for localized fluctuations in smaller 7]-(j) bins 
by comparing with several kinds of mixed events. We give a summary in section 7. 

2. Multiplicity Distributions 

The WA98 experimental setup consists of large acceptance hadron and photon spec- 
trometers, detectors for charged particle and photon multiplicity measurements and calorime- 
ters for measuring transverse and forward energies. The charged particle hits were counted 
using a circular silicon pad multiplicity detector (SPMD) located 32.8 cm from the target 
covering 2.35 < rj < 3.75. The efficiency of detecting charged particles has been deter- 
mined in a test beam to be 99%. The detector is transparent to high energy photons. 
The photon multiplicity was measured using a preshower Photon Multiplicity Detector 
(PMD) ||19[ placed 21.5 meters downstream from the target and covering the pseudora- 
pidity range 2.9 < rj < 4.2. The cluster of hit pads, having total ADC content above a 
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Figure 2. Pseudorapidity distributions of charged particles and photons in Pb induced 
reactions at 158-A GeV on a Pb target. The charged particle distributions shown in 
the left panel are for 0-10%, 10-20%, 20-40%, 40-60% and 60-80% centrality bins. The 
centrality bins for photons are as labelled in the figures. 



hadron rejection threshold are identified as photon-like (A^^-uko)- The photon counting 
efficiencies for central to peripheral cases vary from 68% to 73% [|1^,0. The purity of 
the photon sample in the two cases varies from 65% to 54%. The centrality of the inter- 
action is determined by the total transverse energy (-Et) measured in the mid rapidity 
calorimeter. The centralities are expressed as fractions of the measured total transverse 
energy. The most central corresponds to the top 5% of the minimum bias cross section. 
Extreme peripheral events in the 80-100% range were not analyzed. 

Fig. m shows the minimum bias charged particle and photon-like distributions within 
the full acceptance of the detectors. The distributions corresponding to the centrality 
cuts of 1%, 2% and 5% of minimum bias cross section are superimposed on the figure. 

The pseudorapidity distribution of charged particles and photons at different centralities 
are shown in the left and right panels of Fig. ^ respectively. The data have been corrected 
for geometry and efficiency factors. For photons, the filled symbols represent the measured 
data, and the open symbols are reflections of the filled symbols at r]c.m.{— 2.92). The solid 
histograms show the corresponding distributions obtained from VENUS event generator 
pi. 



3. Scaling of particle production 

The gross features of particle production in nucleon-nucleon collisions and reactions 
of light nuclei are well described in the framework of wounded nuclear model [^]. In 
this model the transverse energy and particle production in p-fA and A+A reactions are 
calculated by assuming a constant contribution from each participating nucleon. While 
a scaling with the number of collisions arises naturally in a picture of a superposition of 
nucleon-nucleon collisions, with a possible modification by initial state effects, the scaling 
using number of participants is more naturally related to a system with strong final 
state rescattering, where the incoming particles loose their memory and every participant 
contributes a similar amount of energy to particle production. The scaling behavior of 
particle production may therefore carry important information on the reaction dynamics 
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Figure 3. Scaling behavior of charged particles and photons. Pseudorapidity density of 
charged particles at midrapidity and integrated number of photons are plotted as functions 
of number of participating nucleons. The solid lines show power-law fit to the data, which 
yields the value of the exponent, a = 1.07 ± 0.05 for charged particles and 1.12 ± 0.03 for 
photons, respectively. 



The number of participants have been calculated using the VENUS model. Pseudora- 
pidity density of A^^ch at midrapidity as a function of the number of participants for Pb+Pb 
collisions are shown in the left panel of Fig. |^. The data points show stronger than linear 
increase (shown as dotted line). A scaling relation can be obtained by fitting the data 
points using C x Np^^^. The value of a is extracted to be 1.07 ± 0.05 p2[. Similarly on 



the right side we have shown the integrated number of photons as a function of number of 



participating nucleons. A similar power law fit yields the value of a to be 1.12 ± 0.03[jT8 
The results for charged particles and photons are consistent with each other within the 
quoted errors. These results show that there is a clear participant scaling violation com- 
pared to a purely linear dependence. The scaling violation might have consequences on 
many other signals, for example, on the J/ ip production. 

4. Multiplicity Fluctuations 

A lot of theoretical interest has been generated on the subject of event-by-event fluctu- 
ations, primarily motivated by the near perfect Gaussian distributions of several observ- 
ables for a given centrality bin. If the distribution of a quantity X is Gaussian, then one 
deflnes the amount of fluctuation by the following: 

where ax is the variance of the distribution. That is, the fluctuation of the distribution 
is the variance squared normalized to the mean of the distribution under consideration. 

In Fig. |I], the distributions of A^ch and A^^-uke are shown for centrality bins of 0-1%, 
0-2% and 0-5%. These curves are very good Gaussians with flts giving x^/ndf to be 
close to unity. It has been observed that making the centrality bin broader beyond 0-5%, 
(from 0-6% and beyond) the distributions deviate from Gaussians. Using the mean and 
variance of the distributions at different centrality bins we have calculated the amount 
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Figure 4. Multiplicity fluctuations of charged particles and photons for various centralities 
within the full coverage of the PMD and SPMD. The width of the centrality bins as 
percentage of minimum bias cross section increases along the x-axis. 
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Figure 5. Multiplicity fluctuations of charged particles and photons for various centralities 
in a common coverage of the PMD and SPMD. The width of the centrality bins as 
percentage of minimum bias cross section remain same (2%) along the x-axis. 



of fluctuations using eqn. (1). These results are plotted in Fig. ^ for charged particles 
and photons. It is seen that the amount of fluctuation increases by increasing the width 
of the centrality bin. This is obvious from the fact that by increasing the width of the 
centrality selection bin, the inherent statistical fluctuations also increases. 

Thus it is very important to control centrality properly in all fluctuation studies such 
that the multiplicity distributions are good Gaussians. Keeping this aspect in mind we 
have used a different set of centrality selection by taking 2% cross section bins, such as, 
0-2%, 2-4%, 4-6%, 6-8%, etc. The resulting multiplicity distributions are good Gaussians 
in nature with x^/ndf between 1 and 1.5. The amount of fluctuations are calculated 
for these type of centrality bins and plotted in Fig. ||. The fluctuations for charged 
particles increase weakly by going from highest centrality towards peripheral whereas 
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Figure 6. The left panel shown the correlation of A^-^-uke and Nch for central (top 
10%) events. The hatched region is for VENUS+GEANT simulation results. The 
right panel shows the Sz distribution for the experimental data is shown, overlaid with 
VENUS+GEANT results (solid histogram). The other curves are generated by incorpo- 
rating 25% and 60% DCC in every event. 

for photons the fluctuations are almost uniform. The fluctuations in photons are higher 
compared to those of the charged particles. This difference may arise because majority of 
photons are decay products of 7r°. These results have to be put in perspective in terms 
of the contributions from known sources, such as, effect of finite multiplicity, fluctuations 
because of impact parameter of the reaction, effect of rescattering, BE correlations, and 
resonance decays. 

5. Charged-neutral fluctuations 

We now turn to event-by-event fluctuations in charged particle and photon multiplici- 
ties. We make a global event characterization in terms of the photon and charged particle 
distributions over the full available phase space of the detectors. The main motivation is 
to search for single large size DCC domains by looking for events which fall far beyond 
the correlation line of A^^-uke — -^ch distributions. 

The correlation plot of A'^-uke and Nch for top 10% central events is shown in the left 
panel of Fig. A strong correlation is seen between charged and neutral multiplicities. 
The fitted line of the distribution is represented by the Z axis. The closest distance 
(Dz) of the data points to the correlation axis has been calculated numerically with the 
convention that Dz is positive for points below the Z-axis. We have chosen to work with 
the scaled variable Sz=Dz/<Jdz in order to compare relative fluctuations at different 
multiplicities. Sz distribution for the data is shown as filled circles in the right panel of 
Fig. 1^. the solid histogram shows the results from the VENUS events passed through the 
GEANT simulation package of WA98 experiment. These events are termed as (V+G) in 
rest of the manuscript. The DCC events are expected to show up as non-statistical tails 
in the distribution of Sz- This is illustrated by the other curves where one introduces a 
known about of DCC fluctuation. Introduction of DCC makes the curves broader. Since 
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we do not see any such event in our data sample, we are faced with the possibihties that 
single-domain DCCs are very rare, very small, or both. Based on a simple DCC model 
calculation we have put upper limits on the frequency of DCC production as a function 
of its size. More details may be obtained from [^]. We discuss search for small size DCC 
in the next section. 

6. Localized charged-neutral fluctuations 

After studying the charged-neutral fluctuations on the full available phase space, the 
next interest would be to search for fluctuations in localized (?7-0) phase space regions or 
domains [0]. One of the major interest would be the search for the formation of small 
size DCC as it is supposed to produce clusters of coherent pions in localized domains. The 
probability distribution of neutral pion fraction in a DCC domain follows the relation : 



P(/) = l/2^/ where f = N^./N^. (2) 

Thus DCC formation in a given domain would be associated with large correlated event- 
by-event fluctuations in the multiplicities of charged particles and photons as majority of 
charged particles consist of charged pions and majority of the photons originate from vr'' 
decays. 

The present analysis uses data from top 5% central events, which corresponds to a 
total of 85K events. The pseudorapidity region common to both PMD and SPMD is 
selected. The acceptance in terms of transverse momentum (px) extends down to 30 
MeV/c, although no explicit selection in px is applied. The experimental results are 
compared to simulated V+G events and various types of mixed events. 

6.1. A simple DCC model 

The effect of non-statistical DCC-like charged-neutral fluctuations has been studied 
within the framework of a simple model in which the output of the VENUS event generator 
has been modified. In this, the charges of the pions within a localized rj-cj) region from 
VENUS are interchanged pairwise (tt+tt" ^ tt^tt^), according to the DCC probability 
distribution as given in equation (2). For the present study, DCC events have been 
generated in the range of 3.0 < 1] < 4.0 with varying intervals in A(j). An ensemble of 
events (henceforth referred to as nDCC events) were then generated by mixing different 
fractions of DCC-like events with those of normal events, as appropriate for different 
probabilities of occurrence of DCC. After allowing the 7r°s to decay, all the particles were 
then tracked through the full GEANT simulation program of WA98 experiment. These 
events are then analyzed using the same analysis methods as of the data. 

Two different analysis methods have been used for the search of non-statistical fluctu- 
ations. The first one is the method of discrete wavelet transformation (DWT) and the 
second one is the A^^_iike-Ar^i, correlation. 

6.2. Discrete wavelet transformation 



The DWT method p6| has been utilized very successfully in many fields including image 



processing, data compression, turbulence, human vision, radar, and earthquake prediction 



25 1 . The beauty of the DWT technique lies in its ability to analyze a spectrum at different 



resolutions with the ability to pick up any fluctuation present at the right scale. Simulation 
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studies by Huang et. al. |]T3[ and Nandi et. al. 0] have shown that the DWT analysis 



could be a powerful technique for the search of localized DCC. 

While there are several families of wavelet bases distinguished by the number of coeffi- 
cients and the level of iteration, we have used the D-4 basis [^. The analysis has been 
performed with the sample function, 



^7-like + ^ch 



calculated at the highest resolution scale, jmax = 5. The spectrum of /' is the input to 
the DWT analysis. The sample function is then analyzed at different scales j by being 
re-binned into 2^ bins. The output of the DWT consists of a set of wavelet or father 
function coefficients (FFCs) at each scale, from j = l,...,{jmax — !)• The coefficients 
obtained at a given scale, j, are derived from the distribution of the sample function at 
one higher scale, j + 1. The FFCs quantify the deviation of the bin-to-bin fluctuations in 
the sample function at that higher scale relative to the average behavior. The distribution 
of FFCs for normal events is Gaussian in nature. However, the presence of localized non- 
statistical fluctuations makes the distribution broader, with a larger rms deviation of the 



FFC distribution |13,15,16|. Comparing the rms deviations of the FFC distribution of 



data and mixed events one can get an idea about the localized fluctuations. 

6.3. Nj — like and A^ch correlations 

The correlation between A^^.^^e and A^ch has been studied in smaller 0-segments by 
dividing the 0-space into 2, 4, 8 and 16 bins. The correlation plots of A^^_iike and A^^ch 
have been constructed for each bin. This is shown in Fig. ^ A common correlation axis 
{Z) has been obtained by fitting the above distributions with a second order polynomial. 
The closest distance (Dz) of the data points to the correlation axis has been calculated 
numerically with the convention that Dz is positive for points below the Z-axis. The dis- 
tribution of Dz represents the relative fluctuations of A^^_iike and A^ch from the correlation 
axis at any given (p bin. In order to compare these fluctuations at different scales in the 
same level, we work with a scaled variable, Sz = Dz/Tms{Dz), where Tms{Dz) corre- 
sponds to VENUS events. The presence of events with localized fluctuations in A^^.^^e 
and A'ch, at a given A0 bin, is expected to result in a broader distribution of Sz compared 
to those for normal events at that particular bin. 

6.4. Mixed events 

In order to search for the presence of fluctuations in the experimental data, it is neces- 
sary to understand all detector related effects by generating different sets of mixed events, 
keeping speciflc physics goals in mind. Properly constructed mixed events should preserve 
all detector effects while removing correlations. Four sets of mixed events are generated 
to provide equivalent event samples as real events. In each type of mixed event the global 
(bin 1) A^^_iike — A"ch correlation is maintained as in the real event. 

The flrst type of mixed events (Ml) are generated by mixing hits on both the PMD 
and SPMD hits separately, with no two hits taken from the same detector. Hits within a 
detector in the mixed events are not allowed to lie within the two track resolution of that 
detector. The second kind of mixed events (M2) are generated by mixing unaltered PMD 
hits of one event with the unaltered SPMD hits of a different event. Two more mixed 
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Figure 7. The scatter plot showing the correlation of A^^-uke and Nch for 1,2,4,8 and 
16 bins in azimuthal angle, </>, for a common coverage of the two detectors, PMD and 
SPMD. The data are for central events, corresponding to top 5% of the minimum bias 
cross section. The scatter plots at different bins are fitted to a curve (shown as Z axis) 
with a second order polynomial. The fit parameters are shown. 



Table 1 

Types of mixed events and how they are used for different physics. 





PMD 


SPMD 


Type of fluctuation to probe: 


Ml 


Mix hits 


Mix hits 


Correlated + Individual 


M2 


Unaltered hits 


Unaltered hits 


Correlated 


M3-7 


Unaltered hits 


Mix hits 


A^^ only 


M3-ch 


Mix hits 


Unaltered hits 


Nch only 



event types are possible which are intermediate between the Ml and M2 mixed events: 
(i) the hits within the PMD are unaltered while the SPMD hits are mixed, this is called 
M3-7, and (ii) the SPMD hits are unaltered while the PMD hits are mixed, this is called 
M3-ch. In table 1 we summarise the construction and usefulness of each of the mixed 
events. 

The behavior of the mixed events may be understood from the sample of nDCC events. 
We have constructed all four different kinds of mixed events from the nDCC events and 
analyzed using both the methods. Here we give the results of the DWT analysis. The rms 
deviations of the FFCs for nDCC events and the different kinds of mixed events produced 
from the nDCC events are shown in Fig. In the absence of DCC-like fluctuation the 
rms values of the various types of mixed events are very close to each other. But the 
V+G rms values come out to be lower than those of the mixed events. This is due to the 
presence of additional correlations between A^^ch and the charged particle contamination 
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0.25 




Probability of DCC events 

Figure 8. The rms deviations of the FFC distributions at j = 1 for simulated nDCC 
events with extent A(j)]:,cc = 90° and for various mixed events constructed from those 
events, as a function of the fraction of DCC-hke events present in the nDCC sample. 



in the A^-Y-uke sample. 

The rms deviations for the Ml-type of mixed events are found to be almost independent 
of probability of DCC-like events. Thus it provides a baseline for studying non-statistical 
fluctuations. The rms deviations of the M2 type of mixed events increase similarly, but 
more weakly, than those of the nDCC events. The rms deviations for the M3 sets of events 
are found to lie between M2 and Ml. Thus, the sequence of the mixed events relative to 
the simulated nDCC events (or data) gives a model independent indication of the presence 
and source of non-statistical fluctuations. The simple DCC model used here results in 
an anti-correlation between A^^-uke and Nch due to the "isospin-fiip" procedure used to 
implement the DCC effect. It also results in non-statistical fluctuations in both N^-\i]^e 
and Nch- Thus the M2 mixed events remove only the N^-ni^e^Nch anti-correlation while 
the Ml mixed events are seen to remove all non-statistical fluctuations and correlations. 
The M3 mixed events give intermediate results because they contain only the A^^_iike 
(M3-7) or (M3-ch) non-statistical fluctuations. 

6.5. Results 

The Sz distributions calculated at 4 and 8 bins in angle are shown in the top panel 
of Fig. ^ for data. Ml and V+G events. The experimental data is broader than the 
simulation and Ml events, indicating the presence of additional fluctuations. The FFC 
distributions extracted from the measured ratio are shown in the bottom panel of 

Fig. 1^ for the experimental data, for Ml-type mixed events (from data), and for V+G 
events. The results are shown for scales j =1 and 2, which carry information about 
fluctuations at 90° and 45° in 0, The FFC distributions of the experimental data are 
seen to be broader than both the mixed and V+G events. The result again suggests the 
presence of non-statistical fluctuations. 

The rms deviations of the Sz and FFC distributions as a function of the number of 



bins in azimuth is shown for experimental data, mixed events, and V+G in Fig. |10[ The 
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Figure 9. The Sz and FFC distributions for 4 and 8 divisions in 0. The experimental 
data, Ml and V+G events are shown by sohd circles, solid histograms and dashed his- 
tograms, respectively. The statistics for data and mixed events are the same, whereas the 
distribution for the V+G events is normalized to the number of data events. 



statistical errors on the values are small and lie within the size of the symbols. The error 
bars include both statistical and systematic errors. The systematic errors include effects 
such as uncertainty in the detection efficiencies, gain fluctuations, backgrounds, binning 
variations and fitting procedures. The total systematic error was obtained as the sum in 
quadrature of the individual error contributions. 

Since the mixed events are constructed to maintain the iV^-uke^iVch correlations for 
the full azimuth (bin 1), the rms deviations of data and mixed events for this bin are 
identical. The difference of the rms deviations between data and V+G for this bin is 
the same as shown in the right hand panel of Fig. ^ The effect of correlations due 
to charged particle contamination in A^^.u^g sample mentioned above (see Fig. |^) has 
been removed by rescaling the mixed event results according to the percentage difference 
between the rms deviations of the Sz and FFC distributions of V+G events and those of 
the corresponding V+G mixed events. 

Comparison of the rms deviations of the experimental data and mixed events have been 
made for all the three different cases. 

• (1) Comparison with Ml : From Fig. ^ one notices that the data points are 
higher compared to those of the Ml events for several bins in 0. The differences 
between the data points and Ml mixed events have been quantified by using the 



quantity, a = y (o"fo«)er)^ + ("^npper)^- 2, 4, and 8bins the values of the Sz rms 
deviations of the data are 2.5cr, 3.0cr, 2.4cr larger compared to those of the Ml 
events, respectively. Similarly the rms deviations of the FFC distributions at 4 and 
8 bins for data are 3.7a and 2.8cr larger than those of the Ml events. At 16 and 32 
bins the results for mixed events and data agree within the quoted errors. 
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Figure 10. The root mean square (rms) deviations of the Sz and FFC distributions for 
various divisions in the azimuthal angle. 



This indicates the presence of locahzed non-statistical fluctuations. This result is 
completely model independent as the comparison has been made to mixed events 
generated from data. The observed non-statistical fluctuations may arise because 
of (a) event-by-event correlated iV-y-uke-^ch fluctuations, (b) individual fluctuations 
in A^^-iike and (c) individual fluctuations Nch_- The source of the fluctuation can be 
obtained by comparing the data with the results from M2, M3-7 and M3-ch. 

• (2) Comparison with M2 : From Fig. |1^ it is evident that the rms deviations of 
the M2 events agree with those of the experimental data within errors for all bins 
in Sz and FFC distributions. This implies the absence of event-by-event correlated 
fluctuations in A^-^-uke versus Nc\i. 

If the amount of DCC-like fluctuations in the experimental data were large, then 
the rms deviations for data would have been larger than those of M2 events. Since 
this is not the case, we have extracted upper limits on the probability of DCC-like 
fluctuations at the 90% confidence level by comparing the results from data with 
those obtained from the nDCC events. Using the DCC-model used in this analysis 
we have extracted upper limits on the probability of DCC events to be 10~^ for 
between 45-90° and 3x10-^ for A0 between 90-135°. 

• (3) Comparison with M3-7 and M3-ch : The M3 type mixed events shown in 
Fig. |10] are found to be similar to each other within the quoted errors and lie between 
Ml and M2. This indicates the presence of localized independent fluctuations in 
AL_iike and A^ch- 
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7. Summary 

Multiplicity distributions and fluctuations are important for understanding the evo- 
lution of nuclear systems at high energy collisions. These have been studied from the 
multiplicity measurements of charged particles and photons in 158-/1 GeV Pb+Pb reac- 
tions. Production of charged particles and photons over the full range of centrality could 
be described in terms of participant scaling as Np^l^'^'^^ and Np^^'^^^^^', respectively. This 
indicates violation of the naive wounded nucleon model. Multiplicity fluctuations have 
been studied by varying the centrality of the reaction. It has been shown that the central- 
ity of the reaction has to be properly chosen so that the multiplicity distributions are of 
good Gaussians. We have shown that this criterion is satisfied by making narrower bins 
in Erp corresponding to 2% of minimum bias cross sections. The fluctuations for photons 
are constant {lo — 2.0 ± 0.06) over the centrality range considered and are higher than 
those for charged particles {u = 1.8±0.11). 

Evcnt-by-event charged to neutral fluctuations have been studied in full acceptance of 
the detectors and in localized {r] — 4>) phase space. Using global event characterization, 
no event with large charged-neutral fluctuations have been observed. A mixed event 
analysis is not possible for this global search. The search for locahzed fluctuations have 
been carried out by comparing data with mixed events generated from the data. Full 
understanding of the nature of the mixed events have been achieved by using a simple DCC 
model. Two different analysis methods, (1) A^-y-ukc and A^ch correlation method, and (2) 
a more sophisticated DWT method, have been employed. Both analysis methods provide 
model-independent evidence for non-statistical fluctuations at the 3a level for intervals 
greater than 45°. This is seen to be due to non-statistical fluctuations in both both 
and Nch- No significant event-by- event correlated fluctuations are observed, contrary to 
expectations for a DCC effect. The origin of the observed individual fluctuations are not 
known at present. The interpretation of the results remains an open question. 

With much higher particle multiphcities achieved at higher energies of RHIC and LHC, 
analysis methods using event-by-event fluctuations would certainly be very essential to 
probe the signals of QGP and DCC. Better understanding of data at SPS energies would 
deflnitely help in making flrmer conclusions of any signal from future experiments. 

REFERENCES 

1. Proceedings of Quark Matter '99, Nucl. Phys A 661 (1999). 

2. Proceedings of Quark Matter 2001, Stoney Brook, USA, to be pubhshed in Nucl. 
Phys. A. 

3. E.V. Shuryak and O. Zhirov, Phys. Lett. B89 (1980) 253. 

4. L. van Hove, Phys. Lett. B118 (1982) 138. 

5. M. Stephanov, et al., Phys. Rev. Lett. 81, 4816 (1998). 

6. Gordon Baym and Henning Heiselberg, Phys. Lett. B469, 7 (1999). 

7. M.Asakawa, U.Heinz, and B.Miiller, Phys. Rev. Lett. 85, 2072 (2000). 

8. S. Jeon and V. Koch, Phys. Rev. Lett. 85, 2076 (2000). 

9. A. A. Anselm, M.G. Ryskin, Phys. Lett. B 266, 482 (1991). 

10. J.D. Bjorken, Int. J. Mod. Phys. A 7, 4189 (1992). J.D. Bjorken, K.L. Kowalski, C.C. 
Taylor, "Baked Alaska" , SLAC-PUB-6109, Apr. 1993. 



15 



11. J. -P. Blaizot and A. Krzywicki, Phys. Rev. D 46, 246 (1992). 

12. K. Rajagopal and F. Wilczek, Nucl. Phys. B 399, 395 (1993); Nucl. Phys. B 404, 577 
(1993). 

13. Z. Huang, et ah, Phys. Rev. D 54, 750 (1996). 

14. T.K. Nayak, WA98 Collaboration, Nucl. Phys. A638, 249c (1998). 

15. B.K. Nandi, G.C. Mishra, B. Mohanty, D.P. Mahapatra and T.K. Nayak, Phys. Lett. 
B449 (1999) 109. 

16. B.K. Nandi, T.K. Nayak, B. Mohanty, D.P. Mahapatra and Y.P. Viyogi, Phys. Lett. 
B461 (1999) 142. 

17. J. Randrup, Phys. Rev. C 62, 064905 (2000). 

18. WA98 Collaboration, M.M. Aggarwal et al., Phys. Lett. B 458, 422 (1999). 

19. M.M. Aggarwal et al., Nucl. Instr. and Meth. A 424, 395 (1999). 

20. K. Werner, Phys. Rep. 232, 87 (1993). 

21. A. Bialas, A. Bleszynski and W. Czyz, Nucl. Phys. B 111, 2000 (1976). 

22. WA98 Collaboration, M.M. Aggarwal et al., Eur. Phys. J CIS, 651 (2001). 

23. WA98 Collaboration, M.M. Aggarwal et al., Phys. Lett. B 420, 169 (1998). 

24. WA98 Collaboration, M.M. Aggarwal et al., e-print archive: [nucl-ex/ 0012004] . 

25. A Craps, IEEE Computational Sciences and Engineering, 2, 50 (1995). 

26. Numerical Recipes, Cambridge Univ. Press , 1998. 



